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We discuss the static fluctuations of the d-wave superconducting order parameter A in CuC>2 
planes, due to quasiparticle scattering by charged dopants. The analysis of two-particle anomalous 
Green functions at T — permits to estimate the mean-square fluctuation S 2 = (A 2 ) — (A) 2 , 
averaged in random dopant configurations, to the lowest order in doping level c. Since A is found 
to saturate with growing doping level while S remains to grow, this can explain the collapse of T c at 
overdoping. Also we consider the spatial correlations (A(O)A(R)) for order parameter in different 
points of the plane. 

PACS numbers: 74.20.-z, 74.62. Dh, 74.72.-h 



I. INTRODUCTION 

Recently, a new piece emerged that physicists should 
insert into the puzzle of high-T c superconductivity 
(HTSC) in metal-oxide materials. The advances in STM 
techniques permitted to obtain nanoscopic resolution of 
their electronic structure, revealing a non-uniform pat- 
tern of local density of states and of superconduct- 
ing (SC) d-wave gap. For instance, in underdoped 
I^S^CaCuOg+a; the islands of ~ 3 nm size of almost 
saturated SC order (sharp coherence peaks in the density 
of states) are surrounded by undersaturated areas (low 
peaks); the pattern also remains non-uniform in the op- 
timally doped material but with undersaturated islands 
over the saturated background. 2J This evidence is cor- 
roborated by the results of precision STM experiments 
on HTSC in the mixed state under applied magnetic field 
H <C H C 2, where no common lattice order for Abrikosov 
vortices was found. 0,0 This experimental challenge in- 
voked a massive theoretical response. [5|-|l0j in partic- 
ular, the scenario of two distinct (and competing) SC 
phases was set forth, supposing segregation of 
an undersaturated /?-phase from the saturated a-phase 
and associating it with nucleation of the pseudogap state 
or the so-called d-density wave state. [ll| However, this 
way of reasoning seems to lead to even harder problems 
with microscopic mechanisms for such strong and irreg- 
ular non-uniformity and it is in odds with the conce pt o f 
stripe phases (another piece of the HTSC puzzle El,!!!) 

Perhaps the most realistic approach to the problem is 
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that based on the evident attraction of charge carriers 
in HTSC to each other and to the ionized dopants. 0- 
These two factors are competitive: the SC 
coupling V leads to formation of uniform states of corre- 
lated pairs of charge carriers, while the impurity potential 
Vi leads to formation of localized states of single carri- 
ers and to local reduction of pairing amplitude . [To| . pl| 
Taking into account that the number of ionized dopants 
equals the number of carriers, the " impurity" effect can 
be well competitive with the pairing for all the levels of 
doping. Below we use the method of Green functions to 
estimate the impurity effect on SC properties of doped 
metal-oxides, especially the amplitude and spatial corre- 
lations of fluctuations of the d-wave order parameter. 



II. GREEN FUNCTIONS AND 
SINGLE-PARTICLE PROPERTIES 

We start from the model Hamiltonian: 
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where the Nambu spinors ip^ = (dk,T) a -k 1) include 
fermion operators of normal quaisparticles with the 2D 
dispersion law£k = 2t(2— cos ak x — cos ak y )— /1, the band- 
width W = 8t, and the chemical potential \x; and fj are 
the Pauli matrices. The gap function Ak = ^i^die 2 ^ — 
includes the BCS-shell restriction by the "Debye en- 
ergy" Ed, and the d-wave factor 7k = cos ak x — cos ak y . In 



the long- wave limit: 7k ' 
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hence the observable coherence 
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peak should lie at A Q b s w irpp A, where po = 4/(ttW) 
is the constant density of states of normal metal in this 
limit. The impurity perturbation Vi produces scattering 
of quasiparticlcs at lattice points p randomly distributed 
with the doping concentration c. The true electronic 
spectrum follows from the Fourier transform of two-time 
fermion Green functions (GF's) 



{(a\b)) s = i e^ + ' i0 >*({a(0,6(0)}}^, (2) 
Jo 

(where (. . .) is the quantum statistical average with the 
Hamiltonian ^ and {...,...} is the anticommutator of 
Heisenberg operators), through the spectral formula: 
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Thus, the Nambu matrix Gk — w^Rivk 
diagonal single-particle GF's and its local value G 
N^ 1 J^k enter the two coupled equations 17i] 
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which define the chemical potential fi and the gap pa- 
rameter A in function of the doping level c and of 
the parameters V, Vi. The standard analysis [Tij gives 

Gk = {e — p, — ^k^3 — Akn — Sk)" 1 , where the self- 
energy matrix can be approximated in the self-consistent 
t-matrix form Ek ~ cV(l — GV) -1 with V = ViT3, then 
the integrals in Eq. 0] yield in the approximate analytic 
formulae: 
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In fact, the first formula in Eq. [3] only applies for 
doping levels c well above the metallization threshold 
Cmet ~ £oPo/2, where e = {po[exp(l/Vip ) - is the 
binding energy of impurity localized state. For c < c me t, 
the system is insulating (so non-superconducting) but 
the chemical potential remains positive, vanishing as 
~ c/[po ln(c me t/c)] at c — > 0. The second formula strictly 
applies at c > EdPo, while at c below this value the pre- 
exponential factor should be changed for l^fe^pj'K. A 
more detailed discussion of the averaged single-particle 
spectrum will be given elsewhere. 



III. TWO-PARTICLE PROPERTIES AND 
FLUCTUATIONS OF SC ORDER 



The non-uniform SC order can be described, using 
the operator of local d-wave order on nth lattice site: 
A n = (^)Ek,k'7ke l(k ~ k,)n a k , T a_kM, with the av- 
erage value N~ 1 J2 n (A n ) = A. Its squared dispersion 
S 2 = N- 1 E„(( A „> ~ ( A n) 2 ) is written from Eq. [Das: 
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It is zero in the doped but uniform crystal (at c > 
but Vi = 0), whereas the non-zero contribution to Eq. 
in the lowest order in Vi comes from two successive 
scatterings: k — > k + q, k' — > k' — q, involved into the 
equations of motion for the momentum non-diagonal two- 
particle (anomalous) GF's. In this approximation, we 
have 



2 ^ cV 2 V 2 A 2 v 
2 ^ 

k,k' 



2 2 0(E k - Ek>) 
7k 7k' 



cV 2 V Po 



(7) 



Accordingly to the written in Sec. ^ the dispersion of 
the observable gap is S b s — t^PPqS ~ (ir/2)y / V pa(c — 
c met ) 3 / 2 Vi, and comparing it with Eq. [5] we conclude 
that this dispersion turns comparable with the aver- 
age gap and can bring the SC order to collapse at 
c ~ (e D /Vi) 2 / 3 (Vp )-V 3 exp(^Ve D p 2 /12). 

Finally, we can estimate the spatial correlation be- 
tween the above mentioned fluctuations, alike Ref. [To| . 
from the analysis of the random function /(n) = 
Em c m ex PH m _ n ) 2 /£o]> wnere c m = 1 for dopant 
perturbed sites, m = p, otherwise c m = 0, and £o is 
the BCS coherence length. Then the spatial decay of the 
correlator of fluctuations (/(n)/(n/)) — (/(n))(/(n/)) w 
c(l — c)cxp[— (m — ii) 2 /2£q] defines their correlation 
length £ c w \/2£o- This implies an estimate for the av- 
erage pinning force due to such non-uniform SC struc- 
ture F p ~ poSobsAobs^Q/a) 2 /t; c . Comparison of this and 
other predictions with the available experimental data 
will be the subject of forthcoming studies. 
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